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Abstract 

A quantum mechanical model is used to derive a generalized Landau- 
Lifshitz equation for a magnetic moment, including fluctuations and dissi- 
pation. The model reproduces the Gilbert-Brown form of the equation in 
the classical limit. The magnetic moment is linearly coupled to a reser- 
voir of bosonic degrees of freedom. Use of generalized coherent states 
makes the semiclassical limit more transparent within a path-integral for- 
mulation. A general fluctuation-dissipation theorem is derived. The 
magnitude of the magnetic moment also fluctuates beyond the Gaussian 
approximation. We discuss how the approximate stochastic description 
of the thermal fleld follows from our result. As an example, we go be- 
yond the linear-response method and show how the thermal fluctuations 
become anisotropy-dependent even in the uniaxial case. 



PACS: 75.10.Jm, 75.30.Gw, 76.60.Es 
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1 INTRODUCTION 



The study of thermally induced magnetization reversal was first carried out 
by Brown (^). His approach was to introduce a noise term into the Landau- 
Lifshitz equation of motion, essentially constructing a Langevin-type equation, 
which we call here the Landau-Lifshitz-Gilbert-Brown equation (LLGB). From 
the LLGB equation the Fokker-Plank equation can be derived that describes 
the time evolution of the probability density distribution of the moment ori- 
entations. Solution of this problem was carried out by Brown for the case 
of an axially symmetric potential and later by Coffey et al (H) for non-axially 
symmetric cases. The numerical solution of the Langevin equation was used by 
Lyberatos et al (^), and has since been applied to the study of magnetization re- 
versal by a number of authors, (^; |^. Recently, Wang et al @ have developed 
an approach introducing a tensor form of the damping constant and applied 
this to the calculation of first mean passage time in the case of an elongated 
grain represented as a chain of coupled particles. Many of these calculations 
are motivated by the need to understand high frequency magnetization reversal 
processes in magnetic recording. The process of reading information currently 
involves GMR sensors, the size of which is continually reducing as recording 
densities increase. This led Smith and Arnctt (j^) to the conjecture that noise 
due to magnetization fluctuations in the read head would be a limiting factor on 
the device size. This is clearly an important problem, which has been further 
developed by Smith (^) and Bertram and co-workers (|^) , who have also studied 
the full micromagnetic description of the problem using an approach in which 
the thermal noise is distributed among the normal modes (|lO|). 

Clearly the introduction of thermal fluctuations in the micromagnetic for- 
malism is important both from the point of view of the physics of magnetiza- 
tion processes and also in relation to important practical problems of magnetic 
recording. Central to all models, both analytical and numerical, is the intro- 
duction of a magnetization fluctuation or an effective fleld via the fluctuation- 
dissipation theorem (FDT) ( |ll|) . The FDT has a strong physical justiflcation, as 
discussed in detail by Landau and Lifshitz however it should be stressed 

that it is strictly valid only for small fluctuations about the local minimum. 
A more serious problem with the use of the Langevin equation in the LLGB 
form is the dissipative term itself, which has no microscopic justification. It 
is clearly important to understand the whole problem of fiuctuations and dis- 
sipation within a first priciples quantum mechanical approach, if the limits of 
the current models are to be established and more fundamental theoretical ap- 
proaches derived. The demand for higher density recording media and faster 
switching rates requires the use of structures on the nanometer scale or less. 
Quantum mechanical effects are then bound to become more and more impor- 
tant to consider in these systems. Effects such as magneto-optical interactions 
may even invalidate the simple damping term that is currently used in the 
Landau-Lifshitz equation. This prompted us to investigate whether the LLGB 
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equation can be recovered from a more fundamental treatment rather than the 
ad hoc approach presently used. Hence it seems natural to ask in what limit 
the LLGB equation can be recovered starting from a quantum model. 

In this paper, we make a small step toward that goal. To address the 
above questions, we take a simple quantum model, that of a single particle with 
large spin interacting with a heat bath and an external magnetic field. The 
spin is taken to be large since we are primarily interested in a semiclassical 
representation of the magnetization vector. This simple model is sufficient to 
allow us to study the effects of thermal fluctuations in many different cases, such 
as magnetization switching in a single domain magnetic particle with uniaxial 
anisotropy, an important problem in magnetic recording physics. The bath is 
taken to be of bosonic nature. Nothing else needs to be assumed to enable us 
to include various mechanisms of interaction between the magnetic moment and 
the environment. We calculate the equations of motion of the magnetization 
and that of an associated fluctuating field in the semi-classical limit. Since our 
interest is mainly in the semi-classical limit, coherent states (CS) are the natural 
choice for the representation of the system. These states have the property of 
minimizing the Heisenberg uncertainty relations. Furthermore, the calculations 
are based on expressing the density matrix of the particle-bath system in terms 
of path-integrals as in the Feynman- Vernon formalism |l^; Clearly this 
method allows a consistent treatment of the magnetization and fluctuations from 
the start. If the thermal field is decoupled from the magnetization, the LLGB 
equation will be shown to correspond to a given choice of density of states of 
the reservoir and of its interaction parameters with the magnetic moment. 

This work is able to provide a different angle from which to discuss the 
LLGB equation and its limitation. We also set a basis against which we can 
examine the discrepancy in the recent calculations of the noise spectrum in 
magnetic recording heads (0; |lO|). Therefore, a treatment of the noise problem 
by a self-consistent method, such as the one presented here, can shed some light 
on why this difference exists. We stress that our results are very general for 
the model considered and no linear approximation is assumed. 

The paper is organized as follows. In section two, we introduce a simple 
model Hamiltonian that can describe dissipation and fluctuations. We linearly 
couple a single domain magnetic particle to an external magnetic field and to 
a Bosonic bath with infinite number of degrees of freedom. In section three, 
we show how to calculate the reduced density matrix elements of the magnetic 
moment. The density matrix elements are naturally expressed in terms of path- 
integrals over the configuration space of the moment In section four, we 
derive coupled equations for the magnetization and fluctuations. We show that 
a general fluctuation-dissipation theorem is satisfied. We also demonstrate how 
to recover the LLGB equation by decoupling the thermal fluctuations, taking 
the high temperature limit and constraining the choice of reservoir. In the 
particular case of LLGB, this corresponds to Gaussian fluctuations and constant 
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dissipation. Similar results have been obtained for the case when the magnetic 
moment is replaced by a harmonic oscillator. This is no surprise since in this 
case, the semi-classical approximation corresponds to a particle with large spin. 
In section five, we compare the classical stochastic treatment to this quantum 
treatment. As an example we calculate the fluctuating field for a single domain 
particle with the external field directed along the easy axis as a function of the 
anisotropy constant. Finally in the last section, we summarize our results and 
state results based on our derived equations when generalized to the anisotropic 
case. 



2 DEFINITION OF THE MODEL 



The model we choose is simple but general enough to include many interesting 
physical situations. It is mainly motivated by the recent work of Safonov and 
Bertram (|T6|). They used a two- level impurity system to simulate relaxation 
effects in a single domain grain. They showed that the damping in their model 
is of the Gilbert form. No fluctuations are considered in their calculation. If 
we consider a collection of spins that are independent, then the magnetization 
vector, M, is a simple sum of these coherent spins, 

M = 5Ms — , (1) 

where S is the spin vector and (^/i^/ft is the gyromagnetic ratio. fj,B is the Bohr 
magneton and V is the volume of the system. In the following we set h = 1, 
gUB = ^ and the density y = 1. From now on, we use the words spin and 
magnetic moment interchangeably. 

We take a single spin S (S^ ^ l) and couple it linearly to a set of oscilla- 
tors and to a constant external field H. The former may represent phonons, 
a time-dependent magnetic field or other Bosonic degrees of freedom. No as- 
sumption will be made about the coupling constants or the density of states of 
the reservoir. The Hamiltonian assumes the following form: 

H = -H ■ S + ^ ojkolak + ^ IkolS- + ^ ll^+ak (2) 

where H is a static external magnetic field. S is the magnetic moment of 
a single particle, and are creation and annihilation operators of the 

reservoir. The coupling constants 7^ may be time-dependent, but will be taken 
as independent of time in the final result. The field H is taken along the z-axis, 
the axis of spin quantization. Coupling the z-component of the vector S to 
the reservoir can be easily added, but it will be omitted in this work. This 
Hamiltonian is sufficient to describe all the desired physics. Using the equation 
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of motion for Sz, it is trivial to see that it is not a constant of the motion and 
hence no hnearization is imphed in this modeL 

The operators are in the Heisenberg representation. The spin operator S 
satisfies the usual commutation relations 

[S2,^±]=0 (3) 

with 

S'^^{S+,S^} + Sl (4) 
where the curly brackets are for anticommutation and 

S+ = Sx + iSy (5) 



S- = Sx - iSy (6) 
while the operators of the reservoir satisfy Bose commutation relations, 

[ afe,a+] = Skk' ■ (7) 

Instead of the usual Fock space representation, we use a CS space representation 
for these operators (18; |9[ ^ . 

For a harmonic oscillator with position , momentum pk and frequency tok , 
the CS are defined as eigenfunctions of the annihilation operator Uk = 

afe|*fe) =*fc|*fc> (8) 



with complex eigenvalues, (|l§|)- These states can also be generated from 
the ground state by applying a displacement operator D (zk) which defines a 
one-to-one correspondence between the complex plane and the oscillator states, 

\zk)=D{zk)\Ok) (9) 

and 

D (zk) = exp {zka'^ - z^ak) . (10) 

CS's form an overcomplete basis and satisfy the minimum uncertainty relation. 
Hence they are the most suitable representation for a semi-classical treatment. 
We also adopt the normalization in (uM 



i^kl^'k) = e-^'^^K (11) 
They also satisfy the following relation, the resolution of the identity operator, 

d<^td^k 



2-Ki 



-e-*^*M<I'fc>('i>fc| = 1- (12) 
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The latter relation is essential for a path-integral representation in terms of 
CS's. 

Similarly for the spin states, we use a CS representation ( pT| p2[ ). They are 
defined by analogy to the harmonic oscillator CS's. The spin components in this 
state satisfy a minimum uncertainty relation, i.e., two of the three components 
commute (^3|). As in the harmonic oscillator 'ground' state |0) is required 

from which to generate all the other states. In this case the state with the 
largest Sz component is taken as the reference state. If the z-axis is taken as 
the quantization axis and if we take = j {j + 1), then by definition, we have 



|0> ^ (13) 

and 

Sz\o)=m (14) 

i.e., the state with the minimum fluctuations ( |l7|). The spin CS's are a general- 
ization of the Holstein-Primakoff construction (^4|) . They are defined in terms 
of deviations from the maximum positive z-component of the spin S 

S..\p) = ij-P)\p) ■ (15) 
The CS's are then constructed using 

Im> = ^^exp(/i5_)|0) (16) 

1 + H 



1 ^ f (2j) 



where ^ is a complex number. Since the configuration space of S is the surface 
of a sphere, it will be clearer to have /i parametrize the surface of a sphere 
through a stereographic projection, 

/i = tan (]-e\ e*"^. (17) 
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In this representation, a CS will be represented by a solid angle Jl 



\n) = \e, if) = l^cos -ej exp |tan 1^-9 j e^'^S^ | |0) . (18) 

A useful property for a path integral formulation is that the unit operator has 
the familiar decomposition in terms of projection operators on all CS's, 



2j + l 
47r 



J dn \n){n\ = i. (19) 
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In this representation, the overlap of two coherent states represents an area on 
a sphere, the surface of which is the configuration space of the spin S. The 
overlap is 



2j 

(20) 



( fi'in ) = (cos ^6* cos -e' + sin -9 sin -e'e'i'^-'^') 1 
^ ' ' \ 2 2 2 2 J 

and its magnitude is 

\{n'\n)\=(^±^X . (21) 



Since we plan to use a path-integral technique, we need to write the expectation 
values of the Hamiltonian in the coherent representation. These expectation 
values follow in turn from those of the operators Sz, S+ and S-. The following 
expectation values are deduced from Eq.(^8|) and Eq. (|2C|), 

) =j(l-cos0), (22) 



n|S'+|n) = jsin^e^^ 



n|S'_|rj ) = jsin^e" 



-ie 



(23) 
(24) 



( n\S\n ) = jn. (25) 

n is a unit vector with angles {0,(p). For j ^ 1, the off-diagonal terms of 
the spin operator are smaller than the diagonal ones by a factor of about ^/j. 
Hence they are negligible in the classical limit. This limit will be implicit in all 
subsequent calculations of the reduced density matrix elements. 



3 REDUCED DENSITY MATRIX ELEMENTS 
OF THE SPIN PARTICLE 

In the following, we make use of CS for both the bath degrees of freedom and 
the magnetic moment. The procedure we follow is by now mostly standard. 
Reference (|2^) (and references therein) provides a general overview of these 
methods and hence we omit most of the intermediate steps in our calculation. 

Bosonic CS's were first used by Langer ( p6| ) to study dissipation and fluc- 
tuations in a superfluid, many-body problem. Starting from the equation of 
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motion of the density matrix of the whole system, p, a Landau- Ginzburg equa- 
tion was recovered in the equiHbrium case and a Fokker-Planck equation in the 
classical limit, 

*|f-[^,P]- (26) 

In a CS formulation, only diagonal elements of the reduced density matrix are 
needed. Instead of starting from the above equation, we can instead start 
from an integral representation of the density matrix elements. This method is 
well known and is based on the Feynman- Vernon formalism (|l3|). This path- 
integral approach is in real-time as opposed to the imaginary-time approach 
in equilibrium thermodynamics. Hence questions like approach to equilibrium 
can be studied within this approach. This method has seen many different 
applications since the Caldeira-Leggett (CL) work (p^. The CL model was 
successful in showing how to recover the Langevin equation by coupling an 
oscillator to a bath of oscillators. It seems natural then to ask if the LLGB 
equation can be recovered by coupling a spin to a bath of oscillators. This 
important question does not seem to have been addressed in the literature. 
Spin coherent states arc the natural language to answer this question. Hence, 
wc formulate the question in terms of CS and use path-integral techniques to 
write the density matrix elements of the system. Use of path-integrals with 
spin CS is not as straightforward as in the case of bosons (pTj), nevertheless it 
is the method of choice in this particular problem. 

The calculation we present below takes into account the correct boundary 
conditions as emphasized in ((^7|)). However we avoid using the more abstract 
holomorphic representation in favor of a more geometric one, i.e., in terms of 
solid angles. The physical space for the Hamiltonian , Eq. (||), is the product 
of the Hilbert space of the spin particle and that of the harmonic oscillators. 



\'^k) ■ (27) 



Using the expectation values of the different operators in the Hamiltonian, we 
get the expectation value of the Hamiltonian in the coherent representation, 

H[$*,$,S] = -H,jcoseit)+Y,^kn{t)<^k{t) (28) 

k 

+J E ^'^'^li^) 0(t)e-^'^^'^ - J J2 ^k^kit) sin ^(t)e*'^(*) . 

k k 

From now on, we normalize the magnitude of all spin vectors by j. The reduced 
density matrix element of the spin particle, pff , is by definition the density 
matrix element of the whole system averaged over the states, |$fc), of the bath, 

Pff, (t) = {Sf, \p{t)\ Sf) (29) 
= JH^n^^kiSr;^ \pit)\Sf;^). 
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where |S/) and |S/') are two arbitrary CS of the spin. 

For simphcity, from now on we use the following notation for the functional 
measure of the Bosonic degrees of freedom, 

J Hs^feS^fc = J S (**,*). (30) 



The calculation of density matrix elements is easily carried out using a path- 
integral representation. The propagator of the Bosonic part can be written in 
terms of a path-integral (pO|) 



*(0)=*i 

t 

+i I dt 



D (**,*) exp K^$*(t)<i>fc(i) (31) 







J2m(.t)dt^k (t)-H(**,*,s) 



L k 



52) 



Running from an initial time, i = 0, to time t, we use a real-time path 
integral to average over all intermediate states. The density matrix element is 
then expressed as an integral in terms of the initial density matrix element of 
the system, 

Pff'it) = j ^{^\<i>) j ^n, j j ^{<^\,^,) j ^{<^i,^2) (33) 

x(S/,$;t I ni,$i;0)(f2i,$i;0 IpI fla, $2; 0) (^2, $2; | Sp,$,t). 



We make no assumption about the initial state of the spin particle. Hence 
we have to calculate a forward propagator, a backward propagator and the 
density matrix element at the initial time. The system is assumed to be at 
finite temperature. Since the Hamiltonian is quadratic, the integrations are 
easily carried out in the stationary-phase approximation. We show a few steps 
in the calculation of the forward propagator. Similar calculations are also done 
for the other two terms in Eq. ( |3^ ) . The forward propagator is first written as 
a path integral 



(S/,$;t I f2i,<i>i;G) = (S/,<f 



(34) 



r2i 



S)Si / I)($t,<I'i)expi^$J,fc(t)<i>i,fe(i) 



-I I dt' 
/o 



;^$j,fc(t')5t'$i,fc(t')-H($t,$i,Si) 



exp{iSwz [Si]} 
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The last term is a geometrical term, the Wess-Zmiiino term (|2^) (and references 
therein) 



[Si] = ds dTSi{s,T)- 
Jo Jo 



rOS^^OS^ (35) 

\ OS OT ' 



where Si (s, t) is a homotopy map between the side (z, r2i) and the side (z, S/) 
(p9|). This term therefore represents the area enclosed by the trajectory of the 
spin vector, and hence there is a corresponding two-form (|2^ . Using Stoke's 
theorem, it can be written in terms of a path integral. It is also well known 
( p2|) that in this form the 1-form that results is physically the potential of a 
magnetic monopole at the center of a sphere. Since a full discussion of the 
topological nature of this term is outside the scope of this paper, we refer to the 
above literature for further details. The bath degrees of freedom are eliminated 
by a stationary-phase evaluation of the integral. The phase is an extremum for 
states that satisfy 

^dt^i.k [t) = (36) 

and a similar equation for <f>J j.. Summations over k are implicit in what follows. 
We have 

«9r$i,fe(r) = ujk^i,kiT) - 7fc(r)5_ (37) 

idr^lkir) = iOk^Uir) ~ lliT)S+. (38) 
The solutions with the correct boundary conditions are 

$^,fc(r) = ^i^kc-'^"^ + 16-'^"^ f dt'e''^''''jkit')S_if) (39) 

Jo 

< T < t. 

At the endpoints, we then have: 

$i fc(f) = $i,fce-*"'=* -h ie-*"'=* / dt'e''^>'''-fk{t')S-it') (40) 



^0 



These solutions are then put back in Eq. (p3|). Similar expressions follow 
from the calculations of the backward propagator. The density matrix element 
at the initial time is calculated with the assumption that the bath is initially 
at equilibrium with the spin. The bath relaxes much faster than the spin, a 
reasonable approximation in many problems in magnetics. In this case the 
density matrix is separable at the initial time. The bath density matrix, ps, 
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is then known and its matrix elements can be written explicitly in terms of the 
Hamiltonian, Hb, of the bath only, 

($i|Pb(0)|$2) = ^($i|e-'5««|$2) (41) 



1 

r/3 



#1 

S)$exp{ $*(/3)$(/3) 



+i / dr [i$*(T)$(r) - W(r)] } 



with periodic boundary conditions 

*(0) = *2 (42) 
**(/?) = 

We find, after applying a stationary-phase approximation to the integral, the 
expression 

($1 |ps(0)| $2) = exp 1^ $^,,$2,^6-''"'= I . (43) 

After integrating out the degrees of freedom of the bath, we are left with only 

integrals over paths in the spin configuration space. The effective action of the 
spin is now complex, as usual with dissipative systems. The reduced density 
matrix element is now given by 

P//'(i) = J^(^ij^^2{^i\Psm^2) (44) 

X / ' SSi / ' DS2 exp Lh, [ dt' {Si,, it') - S2,z it')) 

+iSwz [Si] - iSwz [S2] } X W (Si, S2) 

where the last term is entirely due to the coupling between the bath and the 
spin particle. It is given by 

In W (Si, S2) = - /o dt' /o dt" exp (-zwfe (f - t"))Yk (*') Ik {t") (45) 

X { e it" - t') S2,+ it') S2,- it") 
+ ^l-eit"-t'))Si.+ it')Si.^^it")} 

+ /„* dt' /J dt" exp i-z^k it' -t")) it') Ik it") 
S2,+ it')Si,. it") 

/o dt' Jo dt" exp i-iujk it' - t")) {t') Ik it") 
X { it') S2,- it") + S2,+ it') Si,_ it")} 
-7^1 dt' /o dt" exp i^u;, it' - t")) 7. it') 7^ it") 
X { 5i,+ it") 5i,_ it') + S2,+ it") S2,- it')} . 
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By taking the limit of an infinite number of oscillators, this latter term 
becomes responsible for the appearance of dissipation in this model. After 
calculating the elements of the reduced density matrix, we can now calculate 
its time evolution and find a Fokker- Plank type equation as was done in the 
original work of CL (p^. We choose rather to take the semiclassical limit of 
this expression and see under what conditions, if any, a LLGB equation can be 
recovered. 



4 THE SEMI-CLASSICAL APPROXIMATION 

In this section, we find the equation of motion of the magnetization by calcu- 
lating the most probable configurational paths. This is done by calculating the 
path in the reduced density matrix element for the spin field that has the largest 
weight. Then we show that these paths are really the semiclassical limit of the 
classical paths averaged over the thermal fluctuations in the LLGB equation. 
We also show that the fluctuation-dissipation theorem is satisfied. It reduces 
to the Brown form only in the high temperature limit and only in the linear 
response approximation. This approximation fails when the system is highly 
anisotropic (p^. 

To facilitate the taking of the classical limit we make the change of variables, 

S(r) =i(Si(r) + S2(r)), (46) 
D (r) = Si (r) - S2 (r) . 

The variable D represents the fluctuating field that is coupled to the spin and is 
due to the inherent irreversibility in the system. In terms of these new variables, 
the weight functional W (Si, S2) becomes 

W[S,D] =cxp|/odt7odi"e-*"''*'e*"^*"7^(t')7fc(i") (47) 
X {-^D+ (f) D_ {t") + e [t" - t') S+ [t') D_ [t") 
-{l-Q{t"-t'))D^ [t')S. {t")) 

- /o d^' /o rfi"e-^-'=*'e-'=*" {D^ {t') D_ {t"))} 

where Q (r) is the unit step function 

e(r) = i + -^7' Te-^^. (48) 

In these new variables, the phase InW clearly shows that the action of the 
reservoir results in an extra linear coupling in S and D. Moreover, we now have 
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a quadratic term involving the variable D. This quadratic term is easily seen 
to be real and negative, assuring convergence of the sum over all configurations 
of D. The linear term describing interaction of the fields S and D is imaginary, 
however. In fact, it is such a term that gives rise to dissipation in the energy 
of the spin S. As we will show below, D is the field that is associated with 
the classical random field in the LLGB equation. If the coupling term in 
InW is set to zero, then integrations over D are equivalent to averaging over 
a Gaussian fluctuating field. As will be seen below, this is the limit in which 
we recover the LLGB equation. In this case the bath does not depend on 
S. This approximation fails when the Hamiltonian is not symmetric under the 
interchange of the dynamical spin components. 

To find the semiclassical result for the reduced density matrix, we again 
resort to a stationary phase approximation to the phase of the path integrals 
in Eq.(Q). First we impose constraints on the spin magnitude by introducing 
two Lagrangian multiplier fields, 771 (t) and 772 (t). These are then coupled to 
the spin fields Si and S2, respectively, as follows, 

6{SI{t)~1) = J S)r/i(T) e'Idrviir}{sl(r)-l)_ ^49) 

A similar expression holds for the other spin vector variable, S2, and 772- These 
constraints are then put back in the expression for the reduced density matrix 
element. The phase of the path integral is now a function of four independent 
fields 771, ?72, S and D. Extremizing the phase of these paths gives the semi- 
classical solution in the large spin limit. After solving for the constrained fields 
in terms of the spin fields Si and S2, we write the remaining two equations in 
S and D only, obtaining a generalized form of the LLGB equation 

'^^ ^^'^ = Sx f H + t(^) + T(^)) + X W , (50) 



dt' V 7 4 

and 



dD {t') 
dt' 



D X (h + t(^)) +S X W + D X t(-°^ (51) 



The vectors T(^), T^^' and W are associated with dissipation, thermal 
fluctuations and magnitude fluctuations, respectively. In Cartesian form, they 
are respectively given in terms of the functions J and Fp by 



T(^) (u) = 



I dt'Q {u - t') {J {u - t') (<') - J* (u - t') S+ (t')} 
- /o" dt'e [u ~ t') {J {u - t') S- (f) + J*{u^ t') S+ {t')} 





(52) 



13 



T(^J (u) = 



I jl dt'{{\j {u - t') + Fp (u - t')) (f) 
+ {u-t') + F*p{u-t'))D+ it')} 

- lo dt' { (" - t') +Ffj{u- t')) (f) 
-(ij* {u-t') + F;iu-t'))D+ it')] 



(53) 



and 



W(u) = 



-i /p* dt'e it' -u){Jiu- t') it') -J*{u- t') D+ it')} 
dt'Q it' - u) { J (u - t') it') + J*iu- t') D+ it')} 



(54) 



These vectors are in general non-local in time and hence include memory effects 
in the equations of motion for S and D . This type of behavior is clearly needed 
when the relaxation time of the reservoir is of the same order as that of the spin 
particle. We will not discuss such a situation here. We are mainly interested 
to recover the constant dissipation case. Even though we called Eq.(pO[) and 
Eq.(j^) generalized LLGB equations, it is not yet clear how a Gilbert damping 
term can arise in these equations. However through a careful choice of the 
density of states of the bath, the coupling constants and the initial conditions, 
such damping form can be recovered as shown below. 

To describe dissipation, we take the continuum limit for the bath states. Then 
the spectral functions J and Fp are given by 



J(r-r') 



duj A [lo) I7 {ijj)\^ exp [—iuj (r — r')] 



(55) 



and 



F0 ir - r') = r da; ^J^''] I7 iu;)f exp [-^u; (r - r')] • (56) 
Jo exp [Pu] - 1 

Ff3 is simply the nonzero temperature counterpart of J. A (w) is the density of 
states of the bath. In fact the function, 

Gir' -r) = y{r' -T)+F0iT' -r), (57) 
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is the inverse of the free propagator of the field D. 



The vectors S and D are orthogonal as follows from the constraint equa- 
tions, Eq.(|49|). Note that when D is set to zero, the density matrix becomes 
diagonal but the equation of motion for S will still have an extra term besides 
the precessional term that is due to the external field H. This extra term 
T*^'^) (u) clearly always has a damping effect. We conclude that it is the vector 
S that must be identified with the classical magnetization and that D is the part 
that gives rise to the thermal fluctuations in S. Finally, we observe that the 
last term in Eq.(pO|) can not be recovered in the classical limit. This quantum 
mechanical term is not present in the LLGB equation and is beyond a classical 
linear-response treatment of the problem of fluctuations, ft is of higher order 
in D and temperature-independent. It is easy to see that this term gives rise 
to fluctuations in the magnitude of the spin. These fluctuations can not be 
accounted for classically since the magnitude of the magnetization is assumed 
to be constant. 

One important thing to note from Eq.(^0[) is that the vector T^^-' is complex 
and hence, if fluctuations are present, the equation of motion for S becomes 
complex. The physical interpretation of this equation then becomes obscure at 
this level and may not be used as it stands to get the effective magnetization 
of the particle. Having a complex equation for the extremum path of the 
spin is however expected given that the same result happens in the case of the 
harmonic oscillator (^5|). A solution for the fluctuating magnetization is then 
sought through a direct calculation of the propagators in Eq. (|33|) . 

Now we show that a generalized fluctuation-dissipation theorem is satisfied as 
expected for this system since we started from a closed system and integrated out 
a large part of its degrees of freedom. We will also show that it is the vector D 
that should be regarded as the quantum source of the thermal fluctuations in the 
spin system as treated in LLGB. The vector S is then the physically measurable 
magnetization. To better understand the physical meaning of the field D and 
to recover the standard stochastic description of the thermal fluctuations, we 
introduce yet another field, ^ (t). In Eq. (]33|) , we will replace the l.h.s. of the 
following expression with the r.h.s. 



exp dt' dt"'D+ {t') g [t' - t") D_ {t")j = (58) 

AfJ^^ exp (^-i 1^ dt' 1^ dt"^i (f) ^g-^ [f - t") {t") - z ^* dt'^i it') D, (t')) 

where A/" is a normalization constant. The quadratic term in D is then assumed 
to be a result of an averaging over all configurations of the field Hence the 
path- integral for the reduced density, Eq.(^), is now in terms of three fields. 
The field ^ will now result in a third equation of motion. However to recover a 
thermal field similar to that introduced by Brown (flj), we proceed by assuming 
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that the field ^ is classical, i.e., we ignore its equation of motion. At this point 
the fields ^ and H are treated as non-dynamical fields. Next we solve for the 
magnetization S for a given ^ and only then do we average over all configurations 
of ^ with the quadratic weight that we originally ignored in the solution. In fact 
^ becomes the Brown stochastic field if we take the classical limit, that of high 
temperature. Given these observations, we can now assume that the effective 
thermal field with which the spin is interacting is really nothing more than the 
abstract D field that has been introduced in this calculation of the reduced 
density matrix element. In fact ^ has the following correlation functions 

- 2Si,giT-T') (59) 

1 e Z""" , , fl3^\ 7rA(a;)|7(w)|^ , . , 

= -6w / dujuj coth — ^ ^ " ^ ^' exp [-luj (r - r')] . 

TT Jo V 2 / w 



To recover the correlations of the thermal field assumed in the LLGB equa- 
tion, we simply take the high temperature limit and require that the bath sat- 
isfies the condition, 

= a, (60) 

CO 

where a is a constant. This condition provides the simplest relation between 
fluctuations and dissipation. In this case, the correlation functions for the 
random fleld become simply 

(r) iv (r')) = 25u'akT5 (r' - r) . (61) 

A similar condition arises if we replace the spin degrees of freedom by those 
of a harmonic oscillator ([l^). However at high temperature, as we noticed 
earlier, a large spin can be approximated well by an oscillator. This condition 
is, however, still true even if the bosonic degrees of freedom of the bath are 
replaced by fcrmionic degrees of freedom. 

Finally we consider recovering the LLC equation with the Gilbert form of 
damping. Equations (^0|) and ( ^l| ) are very general as they stand and it is not 
clear if the dissipation has the Gilbert form. To deduce the very special case 
of constant damping with the Gilbert form, we set the fluctuations to zero and 
take the following form for the spectral function J, 

J (r' - r) = iaj^5 (r' - r) . (62) 

After an integration by parts of the term containing T'^'^^ in the reduced den- 
sity matrix clement, Eq.(^), the equation of motion for S, the magnetization, 
becomes simply 

-^ = S(r)x H + « -^--^z (63) 
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The boundary terms in the integration by parts are easily dealt with by a 
renormalization of the measure of the path-integral. Keeping in mind the model 
used to derive this result, this equation reduces to the LLG form only in the limit 
of small deviations from local equilibrium. It is also important to note that the 
choice we made for J, Eq.(|6^, is compatible with the Gaussian approximation 
for the thermal field. Hence the stochastic LLG equation is compatible with the 
FDT at high temperatures. However in this case, the application of the FDT to 
the LLG equation is not so trivial as for the generalized LLGB equations, Eqs. 
( pO|) , (H^). In this case the correlation functions will depend on the dynamics 
of the system, i.e., the symmetry of the Hamiltonian. 

Before we end this section, we make a final comment about the condition, 
Eq.(|6^), by which we recovered the LLG limit. If we assume constant coupling 
constants for the interaction between the spin and the bath, we find that the 
density of states must be linear. For phonons, the density of states is quadratic 
and hence, based on this assumption, can not be the major source of the damping 
constant a. In fact dissipation due to currents is believed to be much larger 
(pO|). Ferromagnetic compounds, such as FeNi, show a complex density of states 
for the non-localised electrons, hence a condition such as that given in Eq.(|60|) 
is representative of many competing mechanisms. It is only the lower part of 
the spectrum that is important for a constant dissipation. In fact in Eq.(p5|), 
the limit of integration can not be taken to be infinite for a real bath. This in 
turn will introduce a new cut-off parameter in condition Eq.(|60|) which will be 
system dependent. 



5 ANISOTROPIC PARTICLES 



Perpendicular recording requires particles with relatively high anisotropy for 
long-term storage purposes. Hence a relevant question to ask is how does 
anisotropy interact with the thermal field. This requires a treatment beyond 
the linear response approach. We treat this question in this last section. We 
limit ourselves to the simplest case; that of uniaxial anisotropy along the z-axis. 
In this case the Hamiltonian of the particle-bath system becomes 

U = -H^Sz - KSl + ^cjfca^Ofc - ^7fca^S'_ - -flS-ak, (64) 

where K is the anisotropy constant. The external field is taken along the easy 
axis. Similarly to the above calculation, we find that the reduced density matrix 
elements in the presence of anisotropy become 
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Pjf'it) = jDnijDn2{ni\p,io)\n2) (65) 

X / 'dSi / 'sSaexpLi?, / dt' { Si,, (f) - 82,, {t')) 

Jill JSf, I JQ 

+iKj^ dt' {Si, {t') - Si, (<')) +i{Swz [Si] - Swz [S2])| W (Si, S2) . 



Clearly the bath influence on the magnetic moment with and without anisotropy 
is the same as before, but now there is coupling between the fluctuating field 
and the spin field that is anisotropy-dependent. Therefore the random field 
becomes K -dependent beyond the linear-response approximation. We show 
below how to calculate the new correlation functions of the fluctuations beyond 
the Gaussian approximation. We define two new anistropy related vectors, 



K5 = z{2KS,) , (66) 

and 

= z{2KD,). (67) 

The equations of motion for the spin field and the fluctuating field become in 
this case 



dt 

and 



P = Sx (h + K5 + t(^) + t(^)) + X (W + VLd) , (68) 



'^SD. = DxfH + Ks + T^'^)) + S X (Ki5 + W) + D X T^^^ . (69) 
dt' V / 



We observe that the additive terms on the right now become K-dependent. 
Hence the fluctuations of the magnitude of the magnetization are anisotropy- 
dependent. The equation of motion for D shows that anisotropy dependence 
can be in the precessional term and hence can be recovered even in the classical 
limit beyond a linear-response approach. We next show that this indeed the 
case. 

From Eq.(|65|), the phase of the path integral, which we denote by F, is 
identified with the effective action for the spin system. The procedure to find 



correlation functions is standard (31). To this action we add two independent 



external sources Qi , Q2. These sources will then be coupled to the fields Si and 
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S2, respectively, and will be used to generate two-point correlation functions for 
the fields S and D. Since in the following we restrict ourselves to equilibrium 
properties, the sources Qi and Q2 will be taken as equal, but still arbitrary. 
This is equivalent to making the external field H arbitrary. Hence H will be 
used instead to generate the correlation functions. Since we are solely interested 
in how the correlation functions depend on the anisotropy constant K, at the 
end we take the limit H 0. We also restrict the discussion to small deviations 
of the magnetization from the z— axis. 

First, we define a new functional Z of the external field H, 



[HI 



DS/DS/, j I)f2iS)n2(ni|p.(0)|rJ2) 
j SSiSSa exp [ r [ H ] ] . 



(70) 



This functional is defined in such a way that its variations with respect to H 
generate all correlation functions of the fields S and D at equilibrium. The 
average value of the field D is clearly given by 



1 



5Z 



Z 5Ha [t) 



(71) 



H ►O 



The two-point irreducible correlation function is similarly given by 



SH^ (t) 6Hp [t') 



H ►O 



{D„{t)Dp{t'))-{D^{t)){Dp{t')) (72) 

From the equations of motion of the fields S and D, we deduce the full equation 
of motion of the function , 



+e''^^ {2K) SxsGx^ [t - t') {Sp (t')> + - t') (*)> 

+-"^^^.(0)^(^.(t)), 

where e"^^ is the antisymmetric unit tensor. In the following the quantum terms 
are neglected and we look for corrections to due only to the anisotropy term. 
We solve Eq.(|73|) by iteration, starting from the free propagator. We find 
that the off-diagonal terms are now non-zero and depend explicitly on K . For 
the 1-2 element of we find that 

^g~^^it-t') = A^S{t-t') (74) 
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where A is a constant proportional to the magnetization. To obtain this equa- 
tion we also assumed that the average of the vector D is zero. This assumption 
is due to the fact that the average of the classical stochastic field, is also taken 
to be zero. 

We define G^^, the Fourier Transform of G^^, using 

Qu = ^ / e-(*-*')^-^i (u;) . (75) 

Solving in frequency Fourier space, we find that 

~ 1 -iAK/T 

The other terms behave similarly as a function of K/T. Since within the 
above assumptions, the correlation functions of the classical thermal field are the 
inverse of those of D-field, we conclude that thermal fluctuations and anisotropy 
behave oppositely. The less the anisotropy the higher the level of fluctuations. 
Higher order corrections can similarly be calculated in this manner. 



6 CONCLUSION 



Using coherent states and a simple quantum mechanical model for a single large 
spin particle, we have shown that a generalized form of the Landau-Lifshitz 
equation can be recovered in the limit of high temperature. In this case the 
damping constant provides relaxation to the local equilibrium state. We have 
also shown how fluctuations give two different contributions to the magnetiza- 
tion. One contribution is magnitude conserving and the other is not. We 
derived generalized equations for the magnetization that include non-local ef- 
fects in the dissipation term and go beyond the simple linear-response approach. 
An important immediate result of this work is the dependence of fluctuations 
on the anisotropy of the system. The LLGB equation is clearly inadequate in 
this respect. However these deflciencies can be corrected by using the right 
correlation functions for the fluctuations. Changing the damping to a tensor 
quantity in the LLG equation (|l0| ) to account for noise correctly is then not 
needed. More complicated couplings, other than the linear coupling considered 
here, naturally induce a tensor character for the relaxation. Most of these con- 
clusions are true even in the classical limit. Garcia-Palacios ( |32| ) treated similar 
questions at the classical level and our results agree in that limit. Generalizing 
our model to the anisotropic case does not change these conclusions. In the 
linear approximation, the results of Safonov and Bertram (10) are reproduced 
from our results by taking the vector D to be the stochastic field. However 
the results of Smith and Arnett (|^) are found by taking the time derivative of 
the vector D to be the random force. Clearly, as seen from Eq.(|5l|), in this 
latter case the correlation functions of the random force can not be taken to be 
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independent of the dynamics of the system. An isotropic fluctuation is sim- 
ply not accurate enough for an asymmetric Hamiltonian. More details on this 
comparison will be communicated elsewhere. 
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